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1.0  statement  of  the  E>roblem 

Recent  developments  in  electro-optical  and  infrared  technologies  have 
brought  a  new  dimension  in  various  areas  of  civil  and  military  applications. 
Laser  beams  or  broadband  infrared  radiation  has  been  used  in  remote  sensing, 
thermal  imaging,  night  vision  devices,  seeker /guidance ,  optical  communication, 
designator,  rangefinder,  and  directed  energy  weapon  systems.  Optical 
countermeasures  such  as  obscurant  smokes  have  also  been  developed.  Such 
electro-optical  and  infrared  systems  have  to  detect/send  electromagnetic  energy 
from/to  a  target  at  a  distance  through  the  atmosjAiere.  A  major  problem  with 
BO/IR  systems  is  that  performance  of  those  systems  is  greatly  affected  by  both 
natural  and  man-made  atmospheric  conditions. 

When  an  electromagnetic  wave  propagates  through  the  atmosphere,  it 
experiences  both  linear  and  nonlinear  effects  such  as  gaseous  and  aerosol 
absorption  and  scattering,  turbulence,  thermal  blooming,  air  breakdown,  etc. 
Intensity  of  electromagnetic  radiation  is  attenuated  due  to  absorption  and 
scattering  by  atmospheric  gases  and  aerosol  particles.  This  transmission  loss 
due  to  absorption  and  scattering  also  causes  reduction  in  contrast  of  the 
target  against  the  background  for  optical  or  thermal  imaging  systems.  Amount 

of  beam  attenuation  depends  upon  beam  wavelength,  pathlength,  and  atmospheric 
conditions  such  as  atmospheric  gas  concentration,  water  content,  aerosol 
concentration,  and  aerosol  size  distribution.  Therefore,  realistic  performance 
tests  of  EO/IR  systems  and  obscurant  countermeasures  must  include  an  accurate 
evaluation  of  atmospheric  transmittance. 

Atmospheric  gases  absorb  selectively  radiation  of  certain  wavelengths 
depending  on  their  molecular  structures.  McClatchey  et  al.  [1]  made  a 
compilation  of  all  absorption  lines  of  the  major  atmospheric  gases.  Strong 
absorbing  gases  are  water  vapor,  carbon  dioxide  and  ozone.  Other  absorbing 
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geises  cure  nitrous  oxide,  carbon  monoxide  and  methane.  Long-range  infrared 
systems  cure  generally  operated  at  one  of  the  atmospheric  "windows"  at  3-4  um, 
4.5-5  ym  or  8-12  ym  where  treinsmission  loss  due  to  absorption  by  atmospheric 
geises  are  relatively  low.  The  doninant  mechanism  causing  transmission  losses 
in  the  atmospheric  windows  is  extinction  by  particles.  Aerosol  extinction  is 
thus  an  in^ortant  factor  in  deteming  transmission  losses  and  it  depends  on 
particle  size,  refractive  index,  particle  shape  and  beam  wavelength. 
Millimeter  waves  have  been  used  to  avoid  problems  resulting  from  beam 
attenuation  due  to  aerosol  particles;  fog,  haze,  etc.  However,  they  are 
affected  by  attentuation  due  to  larger  size  particles  such  as  rain  droplets  and 
snow  flcdces. 

Ihere  have  been  efforts  to  develc^  models  which  can  calculate  spectral 
transmittances  under  different  meteorological  and  atmospheric  conditions.  A 
ocrpjter  code,  LCWTRAN  [2],  has  been  developed  for  calculating  the 

transmittance  and  background  radiance  of  the  earth's  atmosphere  fran  0.25  to 

-1 

29.5  un  at  moderate  spectral  resolution  of  20  cm  .  The  HITRAN  [3],  which  can 
calculate  the  high-resolution  atmospheric  transmittcince  of  the  laser  beams,  has 
been  developed  by  personnel  at  the  Air  Force  Geophysics  Laboratory.  Sane  other 
tTcUTsmission  models  have  also  been  developed.  The  Electro-Optical  Systans 
Atmospheric  Effects  Library  (BOSAEL)  includes  models  and  computer  codes  that 
aid  the  user  in  determining  the  effect  of  atmospheric  and  man-^nade  obscurants 
on  radiation  transmission  [4] .  The  greatest  uncertainty  in  the  application  of 
ccnpiter  models  to  the  on-site  calculation  of  the  transmittance  over  optical 
paths  is  in  the  wavelength  dependent  attenuation  of  electromagnetic  radiation 
by  aerosol  pcurticles. 

Aerosol  attenuation  becomes  very  important  when  an  electro-optical  or  an 
infraured  system  is  used  under  adverse  atmospheric  conditions  such  as  at  low 
visibility  or  when  obscurant  particles  are  present  in  the  path.  The  capability 


of  predicting  performance  of  such  systans  under  various  atmospheric  conditions 
is  essential  for  the  system  designer  as  well  as  the  end  user  of  such  systems. 
It  is  necessary  to  characterize  the  optical  properties  of  individual  aerosol 
particles  in  order  to  predict  performance  of  optical  and  infrared  systans. 

Particles  present  in  the  earth's  environment  from  the  surface  boundary 
layer  to  orbital  altitudes,  which  affect  performance  of  optical  and  infrared 
systems,  range  from  natural  and  artificial  atmospheric  aerosols  through  cloud 
droplets  and  meteoric  dust.  Hygroscopic  aerosol  such  as  cloud,  fog,  and  haze 
droplets  can  be  considered  as  spherical  in  their  shape.  Much  of  the 
atmospheric  aerosols  which  originate  from  soil,  salt,  or  pollen,  artificially 
obscuTcint  pcirticles,  ice  clouds,  and  interplanetary  or  meteoric  dust  in  the 
upper  atmosphere  are  known  to  be  irregularly  shaped.  Quantitative  analyses  of 
scattering  by  atmospheric  aerosols  are  limited  by  the  fact  that  little 
information  is  available  about  the  complex  index  of  refraction.  Refractive 
index  of  particles  is  a  function  of  particle  conposition  and  beam  wavelength. 
No  exact  solution  for  irregularly  shaped  particles  scattering  is  available  at 
present.  Therefore,  it  can  be  conclujed  that  well-structured  experiments  to 
measure  the  light  scattering  pattern  of  individual  particles  are  a  better  way 
of  characterizing  the  optical  properties  of  airborne  peirticles. 

Four  different  approaches  have  been  used  to  date  to  characterize  the 
optical  properties  of  aerosol  particles:  remote  measurement,  analysis  of 
collected  samples,  inferring  optical  constants  from  theoretical  calculations, 
and  in-situ  measurement.  Multichannel  solar  radicmetry  and  Lidar  are  examples 
of  remote  measurements  of  the  optical  properties  of  aerosols.  Particles 
collected  on  filters  can  often  be  examined  by  electron  or  light  microscopy  to 
obtain  their  size  and  shapes.  Absorption  by  the  collected  particles  can  be 
measured  using  integrating  plate  [5]  or  diffusing  plate  [6]  methods. 


Scattering  and  absorption  cross  sections  can  be  calculated  using  Mie's 
scattering  theory  [7,8]  for  spherical  particles  from  the  particle  size 
distribution  data  measured  with  an  optical  particle  sizing  instrument. 
However,  this  approach  has  limitations  resulting  frcm  measuranent  error  due  to 
uncertainties  in  particle  refractive  index  and  particle  shape.  Various  in-situ 
measurement  methods  for  characterizing  the  optical  properties  of  atmospheric 
aerosols  ha^^e  been  developed.  In-situ  measurement  techniques  can  be  classified 
into  two  types;  those  which  measure  on  a  large  number  of  particles 
simultaneously  and  those  which  measure  on  individual  particles  one  at  a  time. 
The  former  type  includes  transmissometer,  nephelcmeter ,  extincticineter ,  and 
photoacoustic  method.  Measurements  on  ensemble  of  particles  result  in  seme 
kind  of  average  quantities  for  a  mixture  containing  particles  whose  optical 
properties  and  composition  might  be  quite  different.  Since  the  advent  of  high- 
power  collimated  lasers,  many  measurements  on  differential  scattering  of  single 
particles  have  been  published.  They  usually  take  one  of  two  approaches; 
levitate  the  particles  in  the  scattering  volume  using  an  electric  field  or 
radiation  pressure,  or  make  rapid  measurements  with  a  rotating 
aperture/detector  or  an  array  of  fixed  detectors. 

To  completely  describe  the  optical  properties  of  non-spherical  particles, 
one  necessarily  needs  to  quantify  all  elements  of  the  Mueller  matrix  [7] .  This 
matrix  relates  the  intensity,  polarization  and  phase  state  of  light  incident 
upon  a  particle  to  those  of  light  scattered,  transmitted  or  reflected  by  a 
particle.  Simultaneous  measurement  of  all  sixteen  Mueller  matrix  elements  is 
considered  as  one  approach  which  can  give  information  concerning  size,  shape 
and  chemical/biological  properties  of  the  scattering  particles.  Sane 
successful  measurements  of  the  Mueller  matrix  elements  have  been  reported  [9] . 
However,  this  method  needs  further  improvements  in  theory  and  experimental 
techniques  in  order  to  provide  rapid  particle-by-particle  analysis  of  the 
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irregularly  shaped  particles  [10],  Even  if  one  could  measure  simultaneously  in 
a  short  period  of  time  all  sixteen  elements,  the  resulting  data  are  extemely 
difficult  to  apply. 

The  normalized  angular  scattering  pattern,  the  so-called  phase  function, 
is  an  inportant  parameter  to  quantify  in  predicting  the  particles'  effect  on 
performance  of  EO/IR  systens.  The  measured  phase  function  can  also  be 
correlated  to  the  physical/chanical  characteristics  of  aerosol  particles. 
Scattering  pattern  of  a  nonspherical  particle  can  be  no  longer  symnetric  along 
the  azimuthal  directions.  This  necessitates  that  scattering  measurements  must 
be  made  simultaneously  at  different  directions,  ideally  over  the  4Tr  steradians 
sphere.  Refractive  index  of  a  particle  changes  with  the  wavelength  of  the 
incident  beam.  Thus,  the  phase  function  need  be  measured  at  different 
wavelengths,  especially  in  the  atmospheric  window  regions.  There  exists  a  need 
to  develop  a  new  instrument  which  can  provide  a  real-time,  in-si tu,  spectral 
measurements  of  scattering  characteristics  of  individual  atmospheric  particles 
of  different  shape  and  composition. 

Our  proposed  approach  is  to  measure  and  record  the  multi-wavelength, 
angular  distribution,  in  three  different  scattering  planes,  of  the  scattered 
light  frcm  individual  particles  passing  through  the  small  scattering  volume. 
All  measurements  have  to  be  done  in  a  real-time  manner  without  altering  sample 
integrity  or  the  physical/cptical  properties  of  the  particles.  During  Phase  I 
efforts,  a  literature  survey  has  been  conducted,  system  requirements  have  been 
developed,  technical  feasibility  of  the  developed  concept  has  been 
investigated,  and  the  design  parameters  have  been  optimized. 


2.0  Riase  ^  Objectives  &  Statement  of  Work 

2.1  Phase  I  Technical  Objectives 

The  specific  objective  of  Phase  I  work  were  as  follows; 

a.  The  design  parameters  were  to  be  identified  and  optimized 
for  a  new  multi-wavelength,  multichannel  nephelometer  system, 
vrfiich  inludes  the  saitpling  inlet,  the  optical  systan  and  the 
data  acquisition  system,  capable  of  in-situ,  real-time 
measurement  of  differential  scattering  patterns  of  individual 
aerosol  particles. 

b.  The  technical  feasibility  of  the  developed  system  concept 
was  to  be  investigated. 

c.  Using  state-of-the-art  techniques,  a  prototype  system  was  to 
be  designed. 


2*2  Phase  I  Task  Implementation 

The  Phase  I  research  contract,  #F49620-84-C-0101,  delineates  specific 
tasks  and  schedules  as  follows: 

1)  Section  B,  0001  RESEARCH 
OOOlAA  The  contractor  will: 

a.  Conduct  a  literature  survey  pertaining  to  measurement 
of  optical  properties  of  aerosol. 

b.  Investigate  the  principles  and  constraints  pertaining 
to  the  cheiracterization  of  the  optical  properties 
of  aerosols. 

c.  Develop  ^stem  requirements  for  determination  of 
optical  properties. 

d.  Optimize  the  design  parameters  for  the  proposed  multi¬ 
wavelength,  multi-channel  nephelcmeter  system 
including  the  sampling  inlet,  the  optical  system  and 
the  data  acq^jisition  system. 
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2)  Section  B,  0002  REPORTS 

Reports  identified  below  shall  be  prepared  in 
accordance  with  Exhibit  A  and  delivered  in  accordance 
with  Section  F. 

0002AA  Final  Report. 

3)  SECTION  F.l,  Period  of  Performance  and  Delivery 

The  supplies  and  services  shall  be  delivered  or  performed 
during  the  followir^  period: 

OOOlAA  15  September  1984  through  14  March  1985 

0002AA  FINAL  REPORT  DUE  13  May  1985 
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3. 0  Summary  of  Phase  ^  Research  Work 

All  objectives  of  Phase  I  work  have  been  accomplished  successfully.  Our 
literature  survey  showed  a  need  for  better  understainding  of  optical  properties 
of  atmospheric  and  man-made  aerosol  particles.  Upon  completion  of  the  Hiase  I 
research,  we  have  concluded  that  development  of  the  proposed  multi-wavelength, 
multichannel  nephelometer  system  is  technically  feasible  using  state-of-the-art 
techniques.  Successful  development  of  such  a  system  will  provide  valuable 
informations  in  broad  areas  of  potential  applications  such  as  design  and 
performance  einalysis  of  electro-qptical  curvd  infrared  systems,  remote  sensing, 
modeling  atmospheric  transmission  of  optical/ infrared  radiations,  optical  or 
thermal  imaging  through  a  horizontal  or  slant  path,  optical  countermeasures,  CB 
agent  detection,  direct  energy  weapon  system,  optical  ccronmunication,  radiation 
transfer  in  the  earth's  atmosphere,  air  quality  and  visibility  monitoring, 
combustion  analysis,  process  engineering,  flow  cytometric  analysis,  multi-phase 
flows  analysis,  etc. 

A  brief  suttmary  of  each  task  performed  during  Phase  I  work  is  presented  as 
follows: 

3.1  Literature  Survey 

Existing  literature  related  to  measurement  of  optical  properties  of 
aerosols  has  been  reviewed  through  veirious  technical  journals  and 
government  publications.  Upon  ccmpleticffi  of  literature  survey,  it  is 
concluded  that  a  new  instrument,  which  can  provide  a  fast-response,  in- 
situ,  spectral  measurement  of  scattering  characteristics  of  individual 
aerosol  particles,  needs  to  be  developed  for  better  understanding  of 
qptical  properties  of  atmospl^ric  aerosols  and  their  effects  on  the 
performance  of  electro-optical  and  infrared  systems.  It  is  found  that  no 
nephelometer  systan  developed  to  date  has  capability  of  fulfilling  such 
measurement  requirements.  Relatively  little  work  has  been  directed  toward 
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in-situ,  real-time  measurement  of  scattering  characteristics  of  individual 
aerosol  particles.  Most  of  the  experimental  work  on  scattering 
characteristics  of  aerosols  have  used  a  nephelcmeter  in  acquiring  data  for 
ensemble  of  randomly  oriented  particles  instead  of  a  single  particle.  All 
differential  scattering  measurements  to  date  on  single  particles  have  been 
made  in  one  scattering  plane,  at  one  visible  wavelength,  and  at  relatively 
low  measuraiient  speeds. 

Besides  measurement  techniques  of  optical  prqperties  of  aerosols, 
literature  survey  on  other  related  areas  such  as  scattering  theories  for 
particles  of  different  shape  and  coitposition,  aerosol  sampling  techniques 
and  optical  particle  sizing  has  also  been  performed. 

3.2  Principles  and  Constraints  Pertaining  to  the  Chcuracterization  of  the 

Optical  Properties  of  Aerosols 


Absorption  and  scattering  by  aerosols  are  the  main  mechanisms  which 
affect  the  performance  of  electro-optical  and  infraired  systems  under 
natural  and  man-made  atmospheric  conditions.  Theoretical  background  and 
measurement  techniques  pertaining  to  the  characterization  of  the  optical 
properties  of  aerosols  have  been  reviewed. 

Four  different  approaches  have  been  used  to  measure  the  optical 
prc^rties  of  atmospheric  aerosols:  1)  remote  measurement  (solar 
extinction  auid  bistatic  lidar),  2)  collected  samples  (absorption  and 
diffusion),  3)  inferring  optical  constants  from  theoretical  analysis,  euid 
4)  in-situ  measurement  (nephelcmeter,  transmissometer ,  and  photoacoustic 
method ) .  Principles  and  constraints  of  each  of  those  methods  have  been 
investigated . 


•-  .V  •- 


Scattering  characteristics  of  aerosols,  which  depend  on  particle 
size,  incident  beam  wavelength,  refractive  index,  particle  shape  and 
orientation,  are  more  ccnplicated  and  important  quantities  requiring 
attention  in  meiny  areas  of  applications.  Nephelcmeters  which  measure 
directly  the  scattering  cheiracteristics  of  individual  particles  or  an 
ensemble  of  particles  in  the  scattering  volume  can  be  classified  into 
three  types  based  on  their  optical  configuration;  fixed-angle 
nephelcmeters,  integrating  nephelcmeters,  and  polar  nephelcmeters. 
Principles  euvd  constraints  of  each  of  the  three  types  of  nephelcmeters 
have  also  been  investigated. 

A  system  concept  for  a  new  nephelcmeter  system  capable  of  in-si tu, 
real-time,  multi-wavelength  measurement  of  the  scattering  characteristics 
of  individual  airborne  particles  of  different  shape  and  ccrposition  has 
been  developed. 

3.3  Development  of  System  Requirements 

Systan  requirements  for  the  proposed  in-situ,  multi-wavelength, 
multichcinnel  nephelcmeter  system  have  been  determined  as  follows: 


-0.3  to  20  um  in  diameter  for 
equivalent  spherical  particles 

-spherical  or  nonspherical 
particles 


1)  Particle  size  range; 


2)  Particle  shape; 


3)  Measurement  angle;  -  5“ '''175®  in  three  different 

scattering  planes 
5 

4)  Particle  concentration;  -  ^  10  particles/cc 

5)  Particle  retractive  index;  -non-absorbing  or  absorbing 

6)  Particle  sanpling  system;  -minimize  size  biasing  during 

particle  sampling 
-provide  single  particle 
measurement 


7)  Measurement  wavelengths;  -three  laser  wavelengths 

(0.6328,  3.39  and  10.6ym) 


8 )  Measurement  speed ; 

9)  Measuranent  accuracy; 

10)  Scattering  intensity 
measurement; 

11)  Data  representation; 

12)  Power  requiranent; 

13)  Environmental  requirement; 


-  ^  100  particles/sec 

-better  than  +  10%  in  forward 
directions 


-two  polarization  ocnponents 
perpendicular  and  parallel  to 
the  scattering  plane 

-store  and  analyze  the 
measured  data  with  a  corputer 
based  fast  data  acquisition 
system 

-provides  real-time ,  in-situ 
measurement 

-110  VAC,  60Hz 

-laboratory  and  field  operable 


3. 4  Optimization  of  the  Design  Parameters 

Technical  feasibilities  of  the  developed  system  concept  to 
satisfy  the  system  requireronts  have  been  investigated.  Based  on  the 
results  of  technical  feasibility  study,  an  in-situ,  multi-wavelength, 
multichannel  nephelcmeter  system  including  the  sanpling  inlet,  the 
optical  system  and  the  data  acquisition  system  has  been  designed. 


a.  Sanpling  Inlet 

The  particle  sampling  system  shall  sanple  airboime  particles  and 

direct  them  through  the  scattering  volume  of  the  receiving  optics  with 

minimal  size  biasing.  Wedding  PM  inlet  for  a  dichotomous  sampler 

10 

initially  developed  for  monitoring  ambient  particulate  matter  shall  be 

used  as  a  sanpling  inlet  for  the  proposed  nephelcmeter  system. 

Effectiveness  of  the  Wedding  PM  inlet  had  been  well  characterized 
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previously  through  a  series  of  wind  tunnel  tests  [11].  It  has  D  value, 

50 

associated  with  a  50%  effectiveness  value,  of  10.7  vim  and  samples  aerosol 
particles  at  a  flow  rate  of  16.67  jipm  with  minimal  dependence  on  wind 
^)eed  and  wind  direction.  The  inlet  employs  a  unique,  omnidirectional 
cyclone  fractionator.  In  dynamic  tests  with  solid  particles,  the  inlet 
displayed  no  particle-boionce  problem.  A  schematic  of  the  particle 
sanpling  system  is  shown  in  Figure  5  in  Section  4.2. 

b.  Optical  System 

Fecisibilities  of  various  optical  configurations  developed  for  single 
particle  differential  scattering  meeisurements  have  been  investigated. 
Theoretical  calculations  on  the  expected  signal  power  levels  incident  on 
detectors  have  been  performed  to  determine  the  measuranent  limits  of  the 
prc^)osed  nephelcmeter.  Lower  limit  of  measurement  size  range  of  the 
nephelometer  design  was  determined  as  0.3  um  for  visible  wavelength  and 
1.0  pm  for  infrared  wavelengths.  An  optical  system  anploying  ellipsoidal 
reflectors  and  rotating  apertures  has  been  designed  for  the  proposed 
nephelometer.  A  schematic  of  the  optical  system  is  shown  in  Figure  6  in 
Section  4.3.  Three  laser  beams  are  focused  on  the  primary  focal  points 
of  circular  strips  of  ellipsoidal  reflectors.  Scattered  light  beams  are 
reflected  by  the  ellipsoidal  reflectors  and  directed  onto  detectors 
located  at  the  secondary  focal  points  of  the  ellipsoidal  reflectors. 
Rotating  apertures  are  placed  in  front  of  detectors  to  make  a  complete 
scan  of  angular  distribution  of  scattered  light  over  2  degrees.  A  more 
oorplete  discussion  of  the  optical  system  is  presented  in  Section  4.3. 


c.  Data  Acquisition  System 

A  fast  response,  microconputer-based  data  acquisition  systen,  which 
can  acquire  the  measured  differential  scattering  data  and  store  them 
simultaneously  into  a  mass  storage  device,  has  been  designed.  A  schematic 
of  the  data  acquisition  system  is  shov/n  in  Figure  7  in  Section  4.4. 
Simultaneous  saiiple  and  hold  technique  is  employed  to  sample  all  nine 
detector  channels  in  a  very  short  time  interval  of  +  5  nanoseconds  in 
order  to  measure  scattered  intensities  corresponding  to  the  same 
scattering  eingle.  The  data  acquisition  system  based  on  an  LSI-11/23 
computer  is  capable  of  anedog  inputs  acquisition  euid  continuous  analog 
througl^t  to  or  from  disk  at  up  to  250,000  samples/sec.  Measured  data 
are  stored  into  a  hard  disk  in  a  real-time  manner  and  transferred  later  to 
a  60  Mbytes  cartridge  t^)e  drive  for  permanent  back-up  storage.  The  data 
acquisition  system  includes  detector  preamplifiers,  signal  conditioning 
electronics,  a  simultaneous  S  &  H  circuit,  an  A/D  converter,  buffers,  a 
corputer,  disk  drives,  a  video  terminal,  a  printer,  and  a  tape  drive. 
Capabilities  of  the  data  acquisition  system  is  further  discussed  in 
Section  4.4. 


4.0  Results  of  Phase  I  Research 

4.1  Review  of  Light  Scattering  Theories  and  Scattering  Calculations 
1)  Scattering  by  Spherical  Particles 

Ihe  scattering  of  an  electromagnetic  wave  by  an  arbitrarily 
shaped  particle  is  usually  described  using  the  Stokes  vector  [7] : 


Index  o  characterizes  the  incident  radiation,  index  s  is  the 
scattered  radiation,  X  is  the  wavelength,  v^ile  R  is  the  distance 
between  the  scattering  particle  and  detector.  S  are  the 

ij 

elements  of  the  4X4  Mueller  matrix.  I  ,  I  »  U  and  V  are  the 

1  2 

Stoke 's  parameters.  All  the  ccnponents  in  the  Mueller  matrix  are 
functions  of  the  scattering  angle,  the  wavelength  of  the  incident 
beam,  and  the  particle  size,  shape,  orientatioi,  and  complex  index 
of  refraction.  Light  plane  polarized  parallel  to  the  plane  of 
scattering  has  Stokes  vector  (1,0, 0,0);  and  light  plane  polarized 
perpendicular  to  the  plane  of  scattering  has  Stokes  vector 
(0, 1,0,0).  The  element  U  and  V  characterize  elliptical 
polarization.  Synmetry  conditions  can  siitplify  the  Mueller 
matrix.  For  randomly  oriented  particles,  Perrin  [12]  has  shown 
that-  the  number  of  independent  elements  reduces  to  six. 
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For  a  single  spherical  particle,  the  Mueller  Matrix  is 


further  reduced  to 


I 

S  0  0  0 

I 

Is 

11 

lo 

I 

0  S  0  0 

I 

2s 

22 

2o 

u 

.  a2 

0  0  s  s 

U 

s 

33  34 

0 

V 

4ir  R 

0  0  -S  S 

V 

s 

34  33 

o 

Mie  [7]  solved  the  Maxwell  equations  related  to  the  scattering 
•  of  light  by  a  homogeneous  dielectric  spherical  particle  with  the 
aE^jropriate  boundary  conditions.  According  to  the  Mie  theory, 
scattered  light  intensity  for  incident  unpolarized  light  is 


1(6) 


i, (a,m,0)+i,(a,m,0) 

lo  - )  (3) 


v^ere  i  and  i  are  Mie  intensity  functions  (=  |E  1^  and  |E  1^;  E 
12  12  1 

and  E  are  ccxtplex  variables);  a  is  size  parameter  and 

2 

m  is  complex  refractive  index  of  particle. 

The  dimensionless  intensity  functions  i  and  i  are 

1  2 

proportional  to  the  square  of  the  electric  field  vectors 

perpendicular  and  parallel  respectively,  to  the  plane  of  the 

observed  scattered  beam.  The  Mie  intensity  functions  are  related 

to  the  Mueller  elements  as 

2  2 


S  =i=lE|,S  =i=|E| 
11  2  2  22  1  1 
S  =  1/2(E  E  *  +  E  E  *) 

33  2  1  1  2 

S  =  i/2(E  E  *  -  E  E  *) 

34  2  1  1  2 


(4) 


For  spherical  particles,  the  values  of  S  ,  S  ,  S  ,  and  S  can 

11  22  33  34 

be  confuted  using  the  computer  code  developed  by  Dave  [13]  once 
the  particle  diameter  and  refractive  index  are  known. 
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The  phase  function,  which  is  defined  as  the  ratio  of  the 
energy  scatterd  per  unit  solid  angle  in  a  given  direction  to  the 
average  energy  scattered  per  unit  solid  angle  in  all  directions, 
is  an  inportant  parameter  to  describe  the  optical  properties  of 
aerosol  particles.  Scattering  by  particles  that  are  very  small 
v\*ien  ccnpared  with  the  wavelength  (o<<l)  is  known  as  Rayleigh 
scattering.  The  Rayleigh  phase  function  is  given  by  [14] 

P^(e)  =  0.75  (l+cos^  6)  (5) 

R 

The  Mie  phase  function  is  given  by 

I 

Ke)  4ttI(9) 

ftn(0)  =  _  =  (6) 

1 

Air 

where  n  is  the  solid  angle  and  C  is  the  scattering  cross  section 

s 

of  the  particle.  A  plot  of  angular  scattering  by  spheres  with 
m=1.55  is  shown  in  Figure  1  [15],  For  small  size  parameters,  the 
Rayleigh  scattering  patterns  are  obtained;  perpendicularly 
polarized  light  is  scattered  isotropically,  while  light  polarized 
parallel  to  the  scattering  plane  vanishes  at  a  scattering  angle  of 
90®.  The  first  deviations  from  Rayleigh  theory  appear  as  forward- 
backward  asyntnetry,  with  more  light  being  scattered  in  forward 
directions.  As  the  particle  size  is  increased  further,  the 


1 


4ir 


i(e)dn 


I  C 
o  s 


asyntnetry  becomes  more  pronounced  and  the  dominant  forward 
scattering  pattern  is  increased  and  more  light  is  scattered  into 
narrow  forward  directons  for  larger  particles. 


scattering  angle 


Figure  1 

Angular  Scattering  by  Spheres  with  ni=1.55;  the  incident 
light  is  polarized  parallel  (-  -  -)  or  perpendicular  (- 
to  the  scattering  plane  [15]. 


For  very  large  spherical  particles,  the  method  of  gectnetric 
optics  is  used  as  an  alternative  to  the  substantial  quantity  of 
computer  time  cind  memory  required  for  a  computation  using  Mie 
theory.  The  contribution  of  diffraction,  reflection  and 
refraction  must  be  calculated  and  added  to  obtain  the  geometric 
phase  functions  [16] . 

2)  Scattering  by  Mon-Spherical  Particles 

No  rigorous  theory  of  scattering  by  irregularly  shaped 
particles  is  available  now  and  in  all  likelihood  the  situation 
will  remain  unchanged  in  the  near  future  even  with  the  advent  of 
high  speed  ccstputers  with  very  large  memory  capacities.  Seme 
scattering  theories  have  been  developed  for  particles  of 
particular  shapes,  which  have  a  syimetry  about  one  or  mere  axes, 
including  circular  cylinders  [17,18],  spheroids  [29,20],  flat 
disks  [21],  fibrous  particles  [22,23]  and  aggregates  of  particles 
[24,25].  Various  analytical  and  numerical  techniques  have  been 
developed  for  the  problem  of  scattering  by  arbitrarily  shaped 
inhanogeneous  particles  such  as:  method  of  moment  [26], 
perturbation  method  [27],  point  matching  method  [28],  extended 
boundary  condition  method  [29-31],  and  semi-empirical  theory  [32]. 
Various  methods  of  calculating  the  light  scattering  properties  of 
nonspherical  particles  have  been  reviewed  by  Schuerman  [33]  and 
Barber  and  Massoudi  [34].  Each  method  has  its  own  advantages  and 
limitations  and  none  of  the  above  methods  can  yield  accurate 
results  for  scattering  by  arbitrarily  shaped  inhomogeneous 
particles. 


Instead  of  the  theoretical  solutions,  empirical  approaches, 
which  are  to  model  the  scattering  occurrences  by  non-spherical 
particles  using  equivalent  shape  (usually  spherical)  as  standard, 
have  been  tried.  Measurements  of  light  scattering  by  numerous 
types  of  irregularly  shaped  particles  [35-38]  indicate  that 
scattering  in  the  near  forvrd  directions  is  similar  to  that  of  an 
ensenble  of  equivalent  spheres.  Microwave  analog  measurements  of 
angular  scattering  patterns  have  been  performed  for  randomly 
oriented  polydisperse  particles  of  different  types:  rough 
spheres,  cubes,  octahedrons,  convex  particles,  concave  particles 
and  fluffy  particles  [39].  Similar  measurements  have  been  done 
for  twenty-eight  particles  of  different  sizes  and  shapes  by 
Schuerman  et  al.  [40].  The  measured  data  were  compared  with  Mie 
calculations  based  on  the  spheres  of  equivalent  volume  or  cross 
sectional  area.  Correspondence  with  prediction  was  good  in  the 
forward  directions  whereas  significant  differences  occurred  at 
medium  scattering  angles  and  towards  backscattering .  Angular 
scattering  patterns  of  irregularly  shaped  particles,  convex  and 
concave  particles,  with  the  size  parameters  in  the  range  of  1.9  to 
5.9  were  similar  to  that  of  spherical  particles  of  equivalent 
volumes.  For  the  larger  particles  whose  size  parameter  varies 


from  5.9  to  17.8,  deviations  in  scattering  intensities  occurred  at 
medium  scattering  angles.  The  intensity  enhancement  observed  for 
backscattering,  which  is  a  typical  behavior  for  mixtures  of 
dielectric  spheres,  is  not  observed  for  nonspherical  particles 
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Pollack  and  Cuzzi  [32]  made  similar  conclusions  for  randomly 
oriented,  irregularly  shaped  particles.  When  their  size  parameter 
is  less  than  5,  their  phase  functions  closely  resomble  those  of 
spheres  of  equal  volume.  For  larger  particles  whose  size 
parameter  is  larger  than  5,  the  phase  function  of  nonspherical 
particles  is  similar  to  that  of  their  spherical  counterparts  at 
small  scattering  angles,  but  shows  marked  deviations  at  larger 
angles. 

Hodkinson  [35]  summarized  the  scattering  behavior  of  a  single 
particle  with  nonspherical  shape.  For  particles  smaller  than  a  = 
0.3,  the  sinple  Rayleigh  scattering  pattern,  whose  absolute 
intensities  are  about  equal  for  particles  of  equal  volume,  is 
observed  regardless  of  their  shape.  With  particles  larger  than 
this,  up  to  ct=  1,  the  forward  lobe  of  the  scattering  pattern  is 
determined  by  particle  dimension  along  the  direction  of  incidence 
and  resembles  the  pattern  of  a  sphere  of  that  diameter.  When 
1  <  a<5  the  scattering  pattern  is  so  cctnplex  that  no 
generalization  with  Mie  calculations  can  be  made.  For  larger 
particles,  when  a  >  5,  the  scattering  can  be  approximated  with 
geometric  optics.  The  angular  spread  of  the  forward  lobe  is 
determined  by  the  transverse  dimensions  of  the  particles,  as  in 
Fraunhofer  diffraction  [41]. 

3)  Scattering  Calculations 

Light  scattering  by  particles  depends  upon  not  only  particle 
size  and  shape  but  also  refractive  index.  Refractive  index  of  a 
particle  is  determined  by  its  chemical  composition  and  is  a 
function  of  wavelength,  too.  Imaginary  part  of  the  refractive 


atnospheric  particles  is  shown  in  Figure  2  [15].  Except  for 
carbon,  there  are  two  distinct  spectral  regions  in  which  k  is  high 
for  those  materials:  one  in  the  infrared  and  the  other-  in  the 
ultraviolet,  with  a  region  of  high  transparency  between.  This 
fact  necessitates  a  need  to  measure  the  optical  properties  of 
atmospheric  aerosols  at  different  wavelengths  of  interest.  Phase 
function  data  measured  at  visible  wavelength  can  not  be  used  to 
predict  and  evaluate  electro-optical  systems'  performance  at 
infrared  wavelengths.  If  particle  size  is  known,  complex 
refractive  index  of  a  particle  can  be  determined  from  angular 
scattering  data  measured  with  a  nephelometer  system.  The  hatched 
region  in  Figure  2  shows  the  approximate  values  of  k  of 
atmospheric  aerosols  at  visible  wavelengths  inferred  by  various 
ranote-sensing  techniques  [42,43].  If  we  compare  these  values  of 
k  with  those  for  individual  constituants  of  the  atmospheric 
aerosols,  it  seems  clear  that  k  determined  remotely  is  some  kind 
of  average  for  a  mixture  containing  a  small  amount  of  a  strongly 
absorbing  component  such  as  carbon  and  much  larger  amounts  of 
weakly  absorbing  components.  The  proposed  nephelometer  systan 
measures  optioal  properties  of  individual  particles  to  get  their 
real  representative  values. 

Total  scattering  cross  section  values  and  angular 
distribution  of  scattered  light  at  different  laser  wavelengths 
have  been  calculated  for  spherical  particles  of  different  size  and 
refractive  index  in  order  to  predict  signal  levels  scattered  into 
different  angles  and  investigate  their  dependency  on  different 
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variables.  A  cxxnputer  code  developed  by  Bohren  and  Huffman  [15] 


was  used  to  calculate  the  scattering  matrix  elements  at  9°  to  180° 
with  a  9°  increment  and  the  scattering  efficient  factor,  Q 

scat 

values.  Values  of  parameters  used  in  scattering  calculations  are: 


wavelength: 


0.442,  0.488,  0.6328,  1.06,  1.15,  3.39,  and 
10.6  ym 


particle  size:  0.3,  0.6,  1,  3,  6,  10,  20,  and  25  ym 

refractive  index:  n  =  1.33,  1.46,  1.55,  1.65  and  1.90 
(m=n-ik)  k  =0.0,  0.1,  0.01,  and  0.005 


Some  results  of  Q  calculations  are  presented  in  Table  1-4. 
scat 


WAVELENGTH  (hw) 


Figure  2 

Imaginary  part  of  the  refractive  index  of  several  solids  and  liquids  that  are 
found  as  atmospheric  particles. 
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Calculation 

for  ^  =  0.6328  Pm 

scat 

0 

Particle 

scat 
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1 - 

refractive  index:  n=1.55- 

- 1 

(pm) 

Parameter 
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0.3 
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10.0 
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1.1513 

20.0 
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2.0576 
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25.0 
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1.1368 
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(pm) 
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0.6 
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1.0 
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2.3266 

3.0 

8.891 

2.8473 

2.5660 

1.2531 

6.0 

17.783 

2.1398 

1.7838 

1.1648 

10.0 

29.638 

2.0366 

1.6196 

1.1593 

20.0 

59.275 

2.0744 

1.4292 

1.1484 

25.0 

74.094 

2.1348 

1.3891 

1.1447 
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Table  3 


Q  Calculation  for  X  =  3.39  ym 
scat 

0 


Particle 

Diameter 

(ym) 

Size 

Parameter 

scat 

n=1.55 - , 

k=U.l 

k=0.0 

“refractive  index: 
k=0.005 

0.3 

0.278 

0.001629 

0.001630 

0.001683 

0.6 

0.556 

0.02644 

0.82643 

0.02704 

1.0 

0.927 

0.1964 

0.1959 

0.1921 

3.0 

2.780 

3.6771 

3.5689 

2.2741 

6.0 

5.560 

2.7247 

2.5692 

1.4617 

10.0 

9.267 

2.8636 

2.5893 

1.2749 

20.0 

18.534 

2.0799 

1.6983 

1.1615 

25.0 

23.168 

2.0906 

1.7554 

1.1635 

Tcible  4 

Q  Calculation  for  X  =  10.6  ym 
scat 

0 


Particle 

Diameter 

(ym) 

Size 

Parameter 

scat 

cc _ 

k=b.O 

“irerractive  inoex: 
k=0.005 

k=0.1 

0.3 

0.089 

0.0001693 

0.00001693 

0.00001751 

0.6 

0.178 

0.0002715 

0.0002715 

0.0002807 

1.0 

0.296 

0.002106 

0.002106 

0.002175 

3.0 

0.889 

0.1682 

0.1678 

0.1656 

6.0 

1.778 

1.6359 

1.6077 

1.2311 

10.0 

2.964 

3.7017 

3.5998 

2.3266 

20.0 

5.928 

2.3376 

2.2100 

1.3150 

25.0 

7.409 

1.9281 

1.7773 

1.1064 
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4.2  Particle 


Sampling  Systan 

The  particle  sampling  systan  must  draw  a  sample  of  atmospheric 
aerosols  with  minimal  size  biasing  and  direct  them  in  a  train  of  single 
particles  into  the  scattering  volume  of  the  systan  optics.  Most  of  the 
optical  instruments  developed  to  date,  including  the  nephelometer  and 
optical  particle  counter,  use  a  sampling  tube  as  a  particle  sampling 
inlet.  This  type  of  sampling  inlet  is  known  to  be  ineffective  because  of 
its  wind  speed  dependency  and  particle  loss  to  the  tube  wall  [44].  To 
obtain  a  representative  measurement  for  particles  present,  one  must  draw  a 
sanple  containing  the  particle  size  range  of  interest  through  an  inlet 
device  into  the  scattering  voliane  of  the  nephelcmeter  system.  To  provide 
meaningful  data,  the  inlet  must  allow  all  particles  of  interest  to  be 
collected  with  the  same  collection  effectiveness  independent  of  sampling 
conditions  in  the  field  environment.  These  conditions  include  mean 
velocity,  turbulence  scale  and  intensity,  and  extraneous  airborne  matter. 

Principal  investigator  and  personnel  of  the  Aerosol  Science 

Laboratory  have  designed,  developed  and  reduced  to  practice  several 

versions  of  Ambient  Aerosol  Inlet  during  the  past  decade.  These  devices 

have  been  exhaustively  laboratory  and  field  tested.  A  review  article  on 

ambient  aerosol  saitpling  has  been  published  by  the  principal  investigator 

[45].  At  the  present  time,  there  are  two  basic  approaches  to  the  sanpling 

of  atmospheric  particulate  matter:  the  Hi-Vol  and  the  dichotomous 

samplers.  The  Wedding  PM  inlet  developed  for  the  dichotorous  sampler 

10 

[11],  which  collects  two  samples  separated  into  fine  and  coarse  (>2.5  ym) 
particle  modes,  will  be  used  as  the  sampling  inlet  for  the  proposed 
nephelcmeter  systen.  The  inlet  is  illustrated  in  Figure  3.  It  employs 
an  omnidirectional  cyclone  fractionator,  allowing  the  aerosol  entry  frcm 
any  angle  of  approach.  An  angular  impetus  is  imparted  to  the  particle 


motion  via  the  eight  evenly  spaced  entrance  vanes.  As  the  particles  enter 
the  inlet  they  follow  the  fluid  stream  lines  along  the  lowsr  radius  and 
enter  the  cyclone  fractionator  through  the  vane  system.  Particle  removal 
is  realized  on  the  inner  collection  tube.  The  flow  with  unremoved 
particles  then  enters  the  middle  tube,  where  the  trajectory  is  altered  and 
the  flow  continues  upward.  An  additional  turn  is  made  to  alter  the  flow 
into  a  downward  trajectory  to  allow  the  particles  to  proceed  to  the 
scattering  volume. 

The  Wedding  PM  Inlet  is  designed  to  be  operated  at  a  sampling 
10 

flow  rate  of  16.67  s,Exn.  Effectiveness  values,  which  is  defined  as  a  ratio 
of  mass  concentration  collected  by  the  inlet  to  that  present  in  the 
ambient  air,  have  been  determined,  as  shown  in  Figure  4,  through  a  series 
of  wind  tunnel  tests  at  three  wind  speeds  of  2,  8,  and  24  km/h.  D 

50 

value,  associated  with  a  50%  effectiveness  value,  of  the  Wedding  PM 

10 

Inlet  is  10.1  '''  10.7  ym.  Effectiveness  value  decreases  to  zero  for 

particles  larger  than  20  ym.  In  order  to  sample  larger  particles, 

modification  of  present  inlet  to  have  D  value  of  15  ym  can  be  considered 

50 

as  an  option  in  the  Phase  II  work. 

A  schematic  of  particle  sampling  system  is  illustrated  in  Figure 
5.  Particles,  after  passing  through  the  inlet,  are  traveling  downward 
inside  the  inlet  exhaust  tube.  In  order  to  achieve  single  particle 
scattering  for  a  high  particle  concentration,  it  is  necessary  to  reduce 
the  particle  velocity  and  to  dilute  the  particle  cloud.  An  isokinetic 
sampling  nozzle  is  placed  at  the  center  of  the  inlet  exhaust  tube  to 
sample  isokinetically  a  part  of  particle  flow  at  a  flow  rate  of  1.0  ipm. 
(Xitlet  flow  from  the  isokinetic  sampling  nozzle  is  expanded  to  reduce  the 
particle  velocity.  A  second  isokinetic  sampling  nozzle  is  used  to  draw  a 
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sample  into  the  scattering  volume  of  the  optical  system.  A  sheath  of 
filtered  clean  air  is  drawn  to  provide  laminar  flow  condition  around  the 
particle  flow  outcoming  from  the  sampling  nozzle.  Clean  sheath  air  is 
also  needed  to  eliminate  the  background  optical  noise  due  to  scattering  by 
particles  in  the  beam  path  outside  the  viewing  volume  of  the  detector 
optics. 

For  optional  aircraft  usage,  a  bag  sampling  system,  which  has  been 
used  successfully  in  particle  size  measurements  with  EAA  (Electrical 
Aerosol  Analyzer)  by  many  workers  [46]  is  proposed.  Real-time  measuranent 
is  difficult  for  aircraft  measurement  because  of  high  flight  speed  of  the 
aircraft.  First  the  vacuum  pump  empties  the  sampling  bag  to  prepare  for 
a  new  sample.  When  the  bag  is  completely  esnpty,  the  operator  closes  the 
bag  evacuation  valve  and  opens  the  aerosol  inlet  valve  until  the  bag  is 
full.  A  sanple  of  air  is  drawn  through  an  isokinetic  sampler  [47]  located 
outside  the  aircraft  boundary  layer.  He  then  closes  the  inlet  valve, 
adjusts  the  flow  rate  through  the  nephelcmeter  systan,  and  begins  the 
measuroTient  cycle. 
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Figure  3  ,  , 

Wedding  PM,ri  Inlet  11 


4.3  Optical  System 


1 )  Background 

Various  optical  configurations  have  been  reported  for  single  particle 
differential  scattering  measuronents.  In  1961,  before  lasers  were  a 
canTon  laboratory  instrument,  Gucker  and  Egan  [48]  published  angular 
light  scattering  measurements  for  single  isolated  particles.  Since  the 
advent  of  high-pover  collimated  lasers  and  high  speed  data  acquisition 
systems,  many  measurements  on  single  particles  have  been  published 
particularly  in  the  last  decade.  To  successfuly  obtain  angular  Light 
scattering  information  from  single  particles  it  is  usually  necessary  to 
take  one  of  two  approaches:  stably  suspend  the  particles  and  take 
measurements,  or  make  rapid  measurements  on  single  particles  in  flow. 
Particle  suspension  methods  can  not  be  used  for  the  proposed  nephelometer 
system  because  of  their  slow  measurement  speed. 

Rapid  light  scattering  measurements,  during  the  short  time  it 
takes  a  particle  to  transit  the  scattering  volume,  can  be  done  in  two 
ways:  with  a  rotating  detector  or  aperture  (or  both),  or  with  an  array  of 
fixed  detectors.  Gucker  et  al.  [49]  used  a  segment  of  ellipsoidal  mirror 
to  direct  light  scattered  by  a  particle  at  one  focus  to  a  detector  fixed 
at  the  other  focus;  an  aperture  {2°  or  5°)  rotating  at  4U00  rpm  then  gives 
a  360®  scan  in  15  msec.  Marshall  et  al.  [50]  used  similar  optical  system 
which  gives  a  360®  scan  in  20  msec.  Morris  et  al.  [51]  mounted  a  detector 
on  a  turntable  that  rotates  at  IHz.  Gram  et  al.  [52]  developd  a  laser 
polar  nephelometer  for  airborne  measurements  which  used  a  detector  rotated 
by  a  stepping  motor  at  a  frequency  of  5Hz.  Moser  [53]  used  two  conical 
reflectors  and  a  rotating  prism  to  scan  360®  in  4  msec.  However,  the  last 
two  systems  measure  light  scattering  by  an  ensemble  of  particles. 
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Fixed  detector  arrays  were  devised  by  Diehl  et  al.  [54],  Crovell  et 
al.  [55]  and  Bartholdi  et  al.  [56,57]  to  measure  angular  scattering  by 
single  particles.  Diehl  et  al.  mounted  detectors  at  +45°,  +9U°,  +  135°, 
combinations  of  either  two  or  three  of  which  sampled  scattered  light  at 
16.7  msec  intervals.  This  discrete  detector  system  can  not  provide  enough 
angular  resolution  in  scattering  measurements.  Crov\ell  et  al.  used  a 
photodiode  array  consisting  of  32  concentric  rings  to  measure  scattered 
light  intensities  simultaneously  at  32  forward  scattering  angles. 
Approximately  250  psec  was  needed  to  aquire  an  event  and  to  transfer  its 
scatter  pattern  to  the  computer.  The  instrument  of  Bartholdi  et  al.  [57] 
has  an  annular  segment  of  an  ellipsoidal  reflector  to  focus  scattered 
light  onto  a  circular  array  of  60  photodiode  detectors.  This  system  was 
designed  for  applications  in  which  biological  cells  flow  in  single  file  at 
a  flow  rate  of  up  to  1000  particles  per  second  through  the  scattering 
volune.  Optical  configuraton  utilizing  a  detector  array  can  provide  high 
speed  scattering  measurements  of  single  particles.  However,  all  systems 
developed  to  date  measure  scattering  intensities  in  one  scattering  plane 
at  one  wavelength. 


2)  Scattered  Radiant  Flux  Calculation 

Expected  radiant  flux  levels  incident  on  the  detectors  have  to  be 

estimated  to  determine  the  signal-to-noise  ratio  and  specific  detector 

types  to  be  used.  Total  radiant  flux  P  ,  scattered  into  4  it  steradian 

s 

sphere  by  a  single  particle  can  be  given  by 

P  =  E  —  ^scat  (7) 

s  \  k 


where  E>  is  irradiance  incident  upon  a  single  particle  in  the  scattering 

2 

volume.  Ex  is  determined  from  R  x  “  Aw  ,  ./here  P  is  light  source 


output  power  and  w  is  a  beam  waist  diameter.  Total  scattered  radiant 

flux,  P  values  have  been  calculated  based  on  0  values  shown  in  Tables 
s  ,  scat 

1-4  for  different  particle  sizes  and  beam  wavelengths,  and  results  are 

shown  in  Table  5.  Refractive  index  value  of  1.55-0.0051,  which  is  known 

to  be  an  average  value  for  atmospheric  aerosols  [58,59],  has  been  used. 

The  beam  waist  diameter,  w,  was  assumed  to  be  100  vim.  Air  molecule  in  the 

scattering  volume  also  scatters  the  incident  light  in  all  directions, 

v^ich  is  Rayleigh  scattering.  The  scattering  volume  size  is  inversely 

proportional  to  sine.  Size  of  the  scattering  volume  averaged  over  4ir 

-5  3 

steradians  was  assumed  to  be  10  cm  ;  100  urn  x  100  um  x  1  mm.  Total 

scattered  power  levels  at  three  laser  wavelengths  due  to  Rayleigh 

scatter ''ig  of  air  molecules  in  the  scattering  volume,  which  are  related  to 

the  background  noise  level,  are  also  calculated  in  Table  5.  Because  of 

small  beam  waist  size,  background  scattering  by  air  molecules  is 

negligible  compared  to  scattering  by  particles  of  interest. 

Radiant  flux  values  shown  in  Table  5  are  not  those  actually  incident 

on  each  detector,  which  would  receive  only  a  small  fraction,  which  is 

subtended  to  the  detector  aperture  opening,  of  total  scattered  light  into 

4tt  sphere.  Angular  distribution  of  scattered  light  also  varies  with 

particle  size  (size  parameter)  as  shown  in  Figure  1.  Power  levels  which 

would  be  incident  on  detectors  located  at  different  scattering  angles  over 

2  TT  degrees  have  been  estimated  and  presented  in  Table  6.  Solid  angle  to 

-3 

which  each  detector  subtends  was  assumed  to  be  6.99  X  10  steradians  for 
Ae=3®  and  &(p=  7.64®.  Total  power  loss  due  to  non-perfect 

transmittances/ref lectances  of  optical  components  throughout  the  beam  path 


From  Table  6,  expected  radiant  flux  values  incident  on  detectors  at 

the  scattering  angles  in  the  forward  to  backward  directions  are  found  as; 

-12  -7  -14  -7  -5  -14 

10  ''' 10  W  at  0.6328  ym,  10  'v  10  W  at  3.39  ym,  and  10  -v  10  W  at 

10.6  um.  Typical  noise  equivalent  power  (NEP)  value  of  a  silicone 
-14 

photodiode  is  10  W/\/hz.  Thus  it  is  feasible  for  a  Si  photodector  to 

measure  the  angular  distribution  of  scattered  light  at  0.6328  ym  by  a 

single  particle  in  the  size  range  of  interest,  0.3 '''20  ym.  However,  in 

the  infrared  wavelengths  no  detector  is  available  for  measurement  of 

-11 

radiant  radiation  less  than  10  watts.  Low  value  of  expected  radiant 

flux  from  a  single  particle  scattering  at  two  infrared  wavelengths 

preclude  the  use  of  uncooled  infrared  detectors.  NEP  of  cooled  detectors 

-12  -11 

is  at  the  order  of  10  '''10  Theormoelectrically  cooled  InAs 

detector  and  liquid  nitrogen  cooled  {77“K)  HgCdTe  detector  are  chosen  for 
3.39  and  10.6  urn  measuranents ,  respectively.  Thus,  at  the  two  infrared 
wavelengths,  3.39  and  10.6  ym,  measurement  of  angular  scattering 
distribution  by  a  single  particle  is  practically  limited  even  with  a 
cooled  detector  to  particles  larger  than  1.0  um.  However,  since  the  size 
parameter  values  of  particles  smaller  than  1.0  ym  are  less  than  one  at 
both  3.39  and  10.6  ym  wavelengths,  scattering  by  those  particles  can  be 
considered  as  Rayleigh  scattering  which  results  in  same  phase  function 
given  by  equation  (5)  independent  of  particle  size,  shape,  and  refractive 
index.  There  is  no  need  to  measure  phase  function  at  3.39  and  10.6  ym 
wavelengths  of  particles  smaller  than  1.0  ym  because  it  is  already  known 
frcm  Rayleigh  scattering. 

Calculations  shown  in  Table  6  were  based  on  scattering  of  a  spherical 
particle  having  refractive  index  of  1.55-0.0051.  Particles  of  different 
shape  and/or  refractive  index  would  result  in  different  angular 
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Table  5 

Total  Radiant  Flux  Scattered  by  a  Single  Particle,  P  (watts) 

s 


Particle 

Diameter 

0. 

,  6328um 

^  -3.39Pin 

^=10.6^^ 

(van) 

(Po=25niw) 

(Po=25inw) 

(Po=3w) 

0.3 

2.01 

X 

-7 

10 

3.67  X 

-10 

10 

4.57  X 

-10 

10 

0.6 

3.24 

X 

-6 

10 

2.38  X 

-8 

10 

2.93  X 

-8 

10 

1.0 

8.69 

X 

-6 

10 

4.90  X 

-7 

10 

6.32  X 

-7 

10 

3.0 

5.12 

X 

-5 

10 

8.03  X 

-5 

10 

4.53  X 

-4 

10 

10.0 

3.99 

X 

-4 

10 

6.47  X 

-4 

10 

1.08  X 

-1 

10 

20.0 

1.27 

X 

-3 

10 

1.70  X 

-3 

10 

2.65  X 

-1 

10 

Air  Molecules 

2.21 

X 

-10 

10 

2.60  X 

-13 

10 

3.27  X 

-13 

10 

Table  6 

Range  of  Expected  Power  Incident  on  a  Detector 

at  0”  '\-180®  Scattering  Angles,  P  (0)  (watts) 

s 


Particle 

Diaineter 


m 

lO 

m 

GIQ 

i3 


X=  0. 6328ym 


p,(e)  — 

A=  3.39viin 


5  X  10 

2  X  10 

-10 

-11 

5  X  10 

10 

-9 

-12 

2  X  10 

6  X  10 

10 

-9 

-12 

10 

10 

2 

X  10 

-7 

-12 

10 

2  X  10 

10 

-7 

-12 

3  X  10 

6  X  10 

4 

X 

O 

A=  lO.Sym 


10  3  X  10 

-5 

2  X  10  1  X  10 

-5 

6  X  10  'V  2  X  10 


distributions/amplitudes  which  may  not  be  much  different  from  those 
presented  in  Table  6.  It  is  concluded  that  it  is  feasible  to  measure  the 
phase  function  of  a  single  particle  at  the  proposed  three  wavelengths; 
0.6328,  3.39,  and  10.6  um.  Three  types  of  detectors.  Si  photodiods,  InAs 

detector,  and  HgCdTe  detector,  were  chosen  for  three  laser  wavelengths, 
respectively. 

3 )  Optical  System 

Three  different  approaches  to  optical  system  design  of  the  proposed 
multi-wavelength,  multichannel  nephelaneter  were  considered  and  their 
feasibility  has  been  investigated  during  Phase  I  effort.  Those  three 
approaches  were: 

I.  use  a  number  of  discrete  detectors  located  over  2 it 
scattering  angles  around  the  scattering  planes. 

II.  use  an  ellipsoidal  reflector  at  the  scattering  plane 
and  detector  arrays  placed  in  front  of  the  second 
focus  of  the  ellipsoidal  reflector. 

III.  same  as  II  but  use  a  rotating  aperture  instead  of 
detector  array  to  scan  entire  scattering  angles. 

Afproach  I  was  initially  proposed  in  Phase  I  proposal.  Optical 
system  is  relatively  sinple  and  easy  to  be  fabricated.  Main  problems  with 
this  approach  are  poor  angular  resolution  due  to  space  limitation  resulted 
from  actual  size  of  detectors,  and  incapability  of  measuring  scattered 
intensities  simultaneously  at  the  same  scattering  angle  at  three 
wavelengths.  It  was  concluded  frcm  scattering  radiant  flux  calculation  in 
the  previous  section  that  three  different  detectors  had  to  be  used  to 
detect  angular  scattered  intensities  by  a  single  particLv,.  A  detector 
such  as  a  pyroelectric  detector,  which  has  a  flat  reponse  over  the  entire 
spectrum  from  visible  to  infrared  regions,  can  not  be  used.  Three 
detectors  have  to  be  placed  closely  at  locations  of  each  scattering  angle 
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to  be  measured.  Actual  size  of  detectors,  especially  large  size  of 
infrared  detectors  including  cooling  devices,  preclude  the  measuring  of 
fine  signature  of  differential  scattering  pattern  for  particles  of 
interest.  In  order  to  measure  scattered  intensities  at  N  points  over  2ir 
degrees,  3N  detectors  cire  needed  for  each  wavelength  to  cover  three 
different  scattering  planes  to  characterize  optical  properties  of 

nonspherical  particles.  In  an  econcmic  point  of  view,  infrared  detectors 
are  expensive  so  that  total  cost  of  the  system  rises  rapidly  as  one 
increases  N.  If  N  is  assumed  to  be  20,  which  is  certainly  not  enough  to 
unfold  fine  shape  of  angular  scattering  patterns  of  a  particle  over  2ir 
degrees,  cost  for  the  total  number  of  detectors  would  be  $390,000. 
Another  drawback  of  this  approach  is  that  it  is  impossible  to  measure  the 
scattered  intensities  at  a  same  scattering  angle  for  three  laser 
wavelengths.  This  makes  it  difficult  to  investigate  dependence  of 
differential  scattering  pattern  on  wavelengths  when  analyzing  data 
measured  with  the  nephelometer.  Thus  the  first  approach  was  withdrawn. 

The  second  approach  is  to  use  an  ellipsoidal  reflector  and  detector 
arrays,  which  is  similar  to  the  syston  employed  by  Bartholdi  et  al.  [57], 
Laser  beams  intersects  the  particle  flow  at  the  primary  focal  point  of  a 
circular  strip  of  an  ellipsoidal  reflector.  The  scattered  lights  are 
reflected  by  the  ellipsoidal  reflector  and  directed  onto  a  circular  array 
of  detectors  placed  in  front  of  the  secondary  focal  point.  Each  detector 
element  senses  signals  emanating  from  a  region  on  the  strip  corresponding 
to  a  range  of  scattering  angles.  The  size  and  positions  of  the  detector 
elements  determine  the  angular  resolution  of  the  differential  scattering 
pattern.  Advantages  of  using  this  approach  are  seen  that  it  is  capable  of 
measuring  the  differential  scattering  patterns  in  very  fine  angular 
resolution,  and  in  a  real-time  manner  with  a  high  speed  data  acquisition 
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system.  Only  one  array  is  needed  for  each  wavelength  per  scattering  plane 
to  scan  over  2ii  scattering  angles.  A  circular  array  of  50  detector 
elements,  v^ose  element  size  is  200  x  200  ym  square,  around  a  0.25" 
diameter  circle  was  considered.  Circular  array  of  photodiodes  for 
measurement  at  visible  wavelength  is  commercially  available.  However,  for 
two  infrared  wavelengths  the  circular  detector  arrays  have  to  be  custom- 
made  ones.  Several  leading  manufacturers  of  infrared  detectors  have  been 
contacted  to  asses  the  feasibility  of  fabricating  such  detector  arrays. 
It  was  found  that  it  is  technically  feasible  to  fabricate  such  circular 
array  detectors  but  cost  would  be  quite  expensive.  A  total  of  six 
infrcired  circular  array  detectors  are  needed  for  two  wavelengths 
measurements  at  three  scattering  planes.  A  very  fast,  sophisticated  data 
acquisition  system  is  needed  to  measure  and  store  the  differential 
scattering  data  signals  frcm  450  detector  channels  in  real-time  mode.  A 
powerful,  high  speed  computer  similar  to  the  HPIOOO  system  has  to  be  used 
to  do  such  a  job.  Signal  conditioning  electonics  including  preamplifier, 
logarithmic  anplifier,  and  peak-sense  electronics  is  also  needed  at  each 
detector  channel  output  to  interface  with  the  input  of  the  data 
acquisition  system.  Total  equipment/parts  cost  of  a  nephelometer 
arploying  such  optical  configuration  was  estimated  to  be  about  $550,000. 
With  this  system  sampling  rate  over  1,000  particles/sec  can  be 
acccmplished.  This  approach  is  technically  feasible  and  can  provide  fast, 
real-time  measursnent  of  differential  scattering  patterns  in  a  fine 
angular  resolution  at  three  wavelengths.  However,  because  of  its 
extremely  high  cost  of  system  equipment,  the  third  approach  was 
considered. 
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The  third  approach  is  to  use  a  rotating  aperture  instead  of  a 
detector  array  to  scan  over  2ir  scattering  angles.  A  schematic  of  optical 
system  of  this  approach  is  shown  in  Figure  6.  A  rotating  aperature  is 
placed  in  front  of  detectors  located  at  the  secondary  focal  point  of  each 
ellipsoidal  reflector.  The  size  of  aperture  limits  the  angular 
resolution,  A0  ,  of  differential  scattering  measurement.  The  aperture  is 
mounted  on  the  shaft  of  a  fast  scanning  motor.  Continuous  scan 
measurements  over  2-rr  degrees  are  possible  with  only  one  detector  for  each 
wavelength  per  one  scattering  plane.  Thus,  a  total  of  nine  detectors  are 
needed  and  this  also  sinplifies  performance  requirements  of  the  data 
acquisition  systan.  In  this  approach,  motor  speed  has  to  be  fast  enough 
to  ccnplete  a  scan  over  2n  degrees  when  a  particle  is  still  moving  inside 
the  scattering  volume.  The  particle  sampling  system  discussed  in  Section 
4.2  was  designed  to  have  a  particle  pass  through  the  scattering  volume  in 
about  3  msec.  A  very  fast  precision  scanning  motor  which  can  be  driven  up 
to  60,000  rpm  is  available.  Ihe  motor  will  be  rotated  at  a  reduced  speed 
of  30,000  rpm,  with  which  it  will  take  2  msec  to  complete  one  scan  over  2tt 
degrees.  Thus,  a  2  msec  long  pulse  emanates  from  each  detector  output 
when  a  particle  passes  through  the  scattering  volume.  A  differential 
scattering  pattern  of  that  particle  can  be  obtained  by  sampling  and 
digitizing  that  pulse  with  a  certain  time  increment.  Total 
acquisition/conversion  time  required  by  the  analog-to-digital  converter  of 
the  data  acquisition  system  is  10  p sec,  which  means  a  maximum  of  200 
samplings  over  2'’^  degrees  is  possible.  Fifty  samplings  over  2tt  degrees 
with  angular  resolution  of  3'*,  which  is  determined  by  the  aperture  size, 
is  proposed.  Sampling  rate  attainable  with  this  approach  is  lower  than 
that  of  detector  array  approach,  but  maximum  sanpling  rate  of  500 
particles/sec  is  still  higher  than  the  system  requirement  which  is  lOU 
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particles/sec.  One  problem  with  this  approach  is  that  the  particle  keeps 
moving  in  the  scattering  volume  while  the  aperture  scans  over  2  it  degrees 
angles.  The  particle's  rotational  motion  in  a  time  period  of  2  msec  can 
be  neglected  under  the  laminar  flow  condition.  Intensity  distribution  of 
focused  laser  beam  across  the  scattering  volume  has  a  gaussian  shape. 
Change  in  scattered  intensity  due  to  uneven  distribution  of  incident  beam 
intensity  across  the  scattering  volume  will  be  corrected  through 
experimental  calibration  of  the  system  with  monodisperse  spherical 
particles  of  known  size  and  refractive  index  and  numerical  analysis  of  the 
measured  data.  Because  of  the  system's  sinplicity  and  econcmical  reasons, 
the  third  approach  anploying  rotating  apertures  has  been  finally  chosen 
for  the  optical  systan  of  the  proposed  multi-wavelength,  multi-channel 
nepheloroeter  system. 

Three  lasers;  a  25mW,  0.6328  ym  HeNe  laser,  a  25n1W,  3.39  ym  HeNe 

infrared  laser,  and  a  3W,  10.6  um  00  laser,  are  chosen  as  light  sources. 

2 

Instead  of  a  1.06  ym  Nd-YAG  laser,  an  infrared  HeNe  laser  emitting  3.39  m 

radiation  is  chosen  because  3.39  ym  falls  in  the  atmospheric  window  region 

of  3  '''  5  y  m  where  many  infrared  systems  are  operated .  Another  reason  is 

the  CW  Nd-YAG  laser  is  bulkier  than  the  HeNe  laser  because  it  needs  a 

cooling  system.  If  the  1.06  ym  laser  is  used,  a  photodiode  can  be  used  as 

a  detector  and  differential  scattering  of  submicron  particles  can  also  be 

measured.  A  sealed  tube  type  CO  laser  is  chosen  for  the  same  reason, 

2 

although  iu  needs  a  cooling  system  and  a  cooled  detector,  because  its 

wavelength  is  in  the  8 12  ym  atmospheric  window  region  and  it  is  widely 

used  in  many  infrared  applications.  A  portable  heat  exchanger  will  be 

used  as  a  cooling  systan  for  the  CO  laser.  Other  lasers  such  as  the  He- 

2 

Cd  laser  (0.442  ym).  Argon  ion  laser  (0.488  or  0.515  ym)  and  turnable 
dye  laser  (0.4'^0.9  ym)  vere  also  considered  as  light  sources.  Use  of 
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three  visible  or  near  infrared  (1.06  ym)  lasers  could  allow  employment  of 
the  circular  array  detector  approach  and  reduce  the  systan  conplexity. 
However,  one  of  the  important  objectives  of  the  proposed  nephelometer 
development  is  to  measure  the  optical  properties  of  atmospheric  aerosols 
over  a  wide  spectrum  range  from  visible  to  far  infrared  wavelengths. 
That  justifies  the  choice  of  two  infrared  lasers  emitting  3.39  and  10.6 
y  m. 

Three  laser  beams  are  focused  with  lenses  onto  the  scattering  volume 
at  the  primary  focal  points  of  ellipsoidal  reflectors.  Anti-reflection 
coated  ZnSe  lenses  are  used  to  focus  two  infrared  laser  beams.  Because  of 
the  refractive  index  change  of  the  optical  material  with  wavelength,  one 
coaxial  beam  approach  consisting  of  three  wavelengths  can  not  be  used. 
Each  of  three  laser  beams  is  focused  with  its  matching  lenses  onto  the 
same  spot  at  the  scattering  volume  using  fine  alignment  adjustments. 

Three  angular  segments  of  ellipsoidal  reflectors  are  mounted  along 
three  scattering  planes  around  a  sphere  centered  at  the  scattering 
volune.  Light  beams  scattered  by  a  particle  at  the  scattering  volume  are 
reflected  by  ellipsoidal  reflectors  and  directed  through  the  rotating 
apertures  onto  detectors  located  at  the  secondary  focal  points.  A 
photosensitive  pickup  is  mounted  across  the  aperture  plate  at  ^ero 
scattering  angle  position.  Its  output,  which  is  used  to  synchronize 
scanning  procedure  of  differential  scattering  measurement,  is  connected  to 
the  data  acquisition  system. 

After  passing  through  the  aperture,  scattered  beams  consisting  of 
three  wavelengths  are  separated  into  three  wavelength  components  using 
ZnSe  hot  mirrors.  In  front  of  each  detector  is  an  interference  filter 
which  prevents  stray  radiation  from  reaching  the  detector  surface. 
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Polarizers  will  be  placed  in  removable  nounts  located  in  front  of  lasers 
and  detectors.  Neutral  density  filters  will  be  inserted  at  the  forward 
scattering  positions  to  reduce  the  beam  intensity  incident  on  the 
detectors.  All  optical  beam  paths  are  enclosed  in  housings  to  keep  the 
background  radiation  level  as  low  as  possible.  Three  types  of  detectors 
will  be  used;  Si  photodiode  for  0.6328  pm,  themioelectronically  cooled 
(195®K)  InAs  detector  for  3.39  pm,  and  a  liquid  nitrogen  cooled  (77“K) 
HgCdTe  detector  for  10.6  pm.  Typical  specifications  of  these  three 
detectors  are  shown  in  Table  7. 


Table  7 

Typical  Specifications  of  Detectors 


Si  Photodiode 

InAs 

HqCdTe 

sensitive  wavelength  (pm) 

0.18'',  1.15 

l.O'v  3.8 

2^  14 

peak  wavelength  (pm) 

0.85 

3.2 

10.6 

responsivity  (A/W) 

0.55 

-14 

1.0 

3  30 

NEP  (W/Vhz) 

-1 

6  X  10 

11 

10 

detectivity,  D*  (onHz  W  ) 

— 

1.0  X  10 

2  X  10 

operating  temperature  (“K) 

ambient 

195 

77 

cooling 

- 

4-stage 

TE  cooler 

liquid  N2 

4.4  Data  Acquisition  System 

A  high  speed,  canputer-based  signal  processing  and  data 
acquisition  system,  which  can  sample  nine  detector  channel  outputs, 
measure  their  amplitudes,  and  store  measured  values  in  the  memory,  is 


needed.  A  2  msec  long  differential  scattering  pulse  frcm  each  of  nine 
detectors  has  to  be  sampled  at  least  50  times  and  a  maximum  of  100 
particles  have  to  be  sanpled  in  one  second.  Thus,  total  sampling  rate  is 
45,000  sanples/sec.  All  these  measured  data  have  to  be  continuously 
stored  on  a  memory  storage  device  at  a  transfer  rate  of  90K  bytes/sec, 
assuming  each  data  consists  of  twD  bytes.  Sanpling  nine  channels  at  25KHz 
per  channel  requires  a  total  aggregate  sampling  rate  of  225  KHz.  A  high- 
perfonnance,  conputer-based  data  acquisition,  which  is  capable  of  analog 
acquisition  to  menory  at  up  to  250,000  samples-per-second  and  continuous 
analog  throughput  to  or  from  disk  at  up  to  100,000  samples-per-second,  has 
been  designed.  A  schematic  of  the  data  acquisition  system  is  shown  in 
Figure  7.  The  system  is  based  on  the  LSI-11/23  microconputer  CPU,  with 
EIS  (Extended  Instruction  Set),  FPU  (Floating-point  Processing  Unit),  MMU 
(Memory  Management  Unit),  256K  Bytes  memory,  and  a  real-time  clock. 

Each  detector  output  is  amplified  with  a  matching  preamplifier  and 
filtered  with  a  bandpass  filter  to  improve  signal  quality.  Logarithmic 
anplifiers  are  used  to  increase  the  dynamic  range  of  measurement  system 
and  improve  system  resolution  for  low-level  signals.  Scattered  radiant 
flux  calculations  in  section  4.1  showed  orders  of  magnitude  change  in 
differential  scattering  data.  Conventional  method  using  a  multiplexer  to 
scan  multi-channel  inputs  is  not  suitable  for  the  present  nephelometer 
application.  Ideally,  one  has  to  sample  nine  channels  at  exactly  the  same 
moment  to  get  scattered  intensity  values  corresponding  to  same  scattering 
angle.  Unavoidable  time  delay  between  each  conversion  when  using  a 


conventional  multiplexer  makes  simultaneous  measurements  impossible.  To 
solve  this  problem,  a  simultaneous  sample  and  holding  (SS&H)  circuit  is 
used  at  the  analog  input  stage  of  the  data  acquisition  system.  This 
method  lets  the  circuit  acquire  nine  analog  input  signals  using  nine 
separate  sample  and  hold  circuits,  which  switch  between  sampling  and 
holding  simultaneously  to  within  +  5ns.  Thus,  the  SS&H  circuit  takes  a 
snapshot  of  nine  analog  input  channels,  freezing  their  values  within  a  +  5 
nanosecond  aperture  uncertainty  interval.  Subsequent  analog  to 
conversions  with  12-bit  resolution  occur  at  up  to  250  KHz  throughput  rate. 

The  gathering  of  data  with  a  microcomputer  based  data  acquisition 
systan  has  been  limited  traditionally  either  to  slow  data  rates  less  than 
10  KHz  aggregate  or  to  small  quantities  of  data  points  fewer  than  16,000 
samples.  Attempts  to  gather  large  numbers  of  samples  at  aggregate  rates 
greater  than  10  KHz  have  resulted  in  data  files  which  contain  gaps  or 
discontinuities  of  an  indeterminate  nature.  In  a  typical  system,  data  is 
output  directly  fron  a  high  performance  DMA  (Direct  Memory  Access)  A/D 
interface  to  memory  and  subsequently  rewritten  from  memory  to  disk. 
Memory  is  typically  organized  as  two  buffers  so  that  while  one  buffer  is 
being  filled,  the  other  buffer  is  being  emptied  to  disk.  Each  word  of 
data  is  tranf erred  over  the  systan  bus  twice,  once  frcm  the  A/D  to  memory 
and  once  from  memory  to  disk.  When  the  data  acquisition  rate  is  greater 
than  10  KHz,  data  are  lost  between  buffers.  A  dual  ported  data 
acquisition  architecture  is  used  in  Figure  7  to  remove  this  constraint. 
The  external  path  connecting  two  boards  allows  the  A/D  board  to  send  data 
to  manory  without  using  the  processor  bus;  the  system  bus  is  needed  only 
to  transfer  data  frcm  memory  to  disk.  The  dual  port  RAM  board  includes  a 
dual  channel.  On-board  Address  Ccxitroller  (OAC)  which  allows  the  memory  to 
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be  configured  as  two  separate  buffers.  The  OAC  can  be  set  up  to  alternate 
between  two  separate  buffers.  The  memory  buffers  are  sufficiently  large 
and  the  buffer  control  logic  is  sufficiently  fast  so  that  this  buffer 
switching  process,  known  as  channel  chaining,  can  be  continued 
indefinitely,  without  loss  of  any  A/D  data  points,  at  rates  up  to  1  MHz, 
subject  to  the  transfer  rate  limits  of  the  A/D  device,  the  processor  bus, 
and  the  mass  storage  subsystem. 

Ihe  data  acquisition  syst«n  includes  two  mass-storage  devices,  a 
floppy /hard  disk  drive  which  provides  1  Mbytes  of  diskette  storage  and  a 
35.6  Mbytes  Winchester  hard-disk  drive.  Hard-disk  drives  are  needed  to 
store  measured  data  at  a  high  transfer  rate  in  a  real-time  manner.  35.6 
Mbytes  of  data  is  equal  to  differential  scattering  data  of  as  many  as 
3,955  particles.  Since  the  disk  in  the  Winchester  drive  can  not  be 
removed,  a  60  Mbytes  cartridge  tape  drive,  which  is  much  slower  than  hard¬ 
disk  drive,  is  used  as  a  back-up  device  for  permanent  storage  of  measured 
data.  Once  the  disk  in  the  Winchester  drive  is  filled  with  data  words, 
measurement  stops  and  the  stored  data  are  dumped  to  the  cartridge  tape  for 
future  data  analysis.  Other  auxiliary  devices  of  the  data  acquisition 
system  are  a  video  terminal  and  a  dot-matrix  printer. 

The  system  software  includes  PT-ll  operating  system,  FORTRAN 
ccrpiler,  and  CPLIB  and  DTLIB  subroutine  packages.  CPLIB  and  DTLIB  support 
analog  input/output,  transfer  of  data,  hardware  interfaces,  and  high 
resolution  plotting.  CPLIB  and  DTLIB  consist  of  a  library  of  subroutines 
designed  to  be  linked  to  user  FORTRAN  programs  and  operated  under  RT-11 
real-time  operating  system. 
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Potential  ccnmercial  applications  of  the  proposed  multi-wavelength, 
multichannel  nepheloroeter  are  seen  for  researchers  and  engineers  in  the 
field  of; 

(1)  Health  Effects  and  Environmental  Instrumentation: 


1.1  Air 


monitoring ;  Inhalable  particulate  matter  whose 


size  is  less  than  15  wm  are  important  for  health  effects 
analysis  because  these  particles  deposit  in  the  tracheo¬ 
bronchial  and  alveolar  regions  of  the  human  respiratory  tract. 
Information  on  the  size  distribution  of  ambient  airborne 
particles  in  that  size  range  is  of  great  importance  in  assessing 
air  quality  and  health  effects. 

1.2  Visibility  degradation ;  Absorption  and  scattering  of 
lightwaves  by  particulates  in  the  atmosphere  is  one  of  the  major 
mechanisms  causing  degradation  of  atmospheric  visibility. 
Characterization  of  optical  properties  of  atmospheric  particles 
with  the  proposed  nephelaneter  system  will  enhance  ones  ability 
of  understanding  effects  of  atmospheric  particles  on  visibility 
degradation. 

(2)  Medical  and  Biological  Research; 

Optical  analysis  of  the  biological  cells  is  a  valuable  tool 
in  medical  and  biological  research.  Light  scattering  technique 
can  be  used  for  the  rapid,  quantitative  analysis  of  individual 
microscopic  objects,  such  as  biological  cells.  Morphology  of 
cells  can  be  classified  from  differential  light-scattering  data 
measured  with  the  proposed  nephelometer  system. 


(3)  Chanical  and  Process  Engineering; 

3.1  On-Line  character i zat ion  of  irregularly  shaped  suspended 
particles  in  multi-t^ase  and  multi-component  system;  Optical 
particle  sizing  instruments  available  at  present  give  erroneous 
results  for  non-spherical  particles.  Measurements  of  an  entire 
differential  scattering  pattern  of  individual  particles  with  the 
proposed  nephelcnveter  will  provide  more  accurate  information  on 
size  and  shape  of  suspended  particles, 

3.2  Canbustion  process  analysis;  In-situ  measurements  of 
particle  loading  and  size  distribution  can  provide  important 
data  for  the  reliable  design  of  canbustion  process  equipment  and 
control  of  environmental  factors.  Examples  of  current  and 
potential  applications  of  the  proposed  nephelcmeter  include  the 
study  of:  ash  fonnation  from  pulverized  coal,  mineral  carry 
over  from  fluidized  bed  combustors  and  pressurized  gasifiers, 
droplet  formation  and  distribution  in  spray  combustors,  and 
radiative  heat  transfer. 

(4)  Radiation  Transfer  and  Atmospheric  Optics; 

The  presence  of  suspended  particulate  matter  in  the  atmosphere 
modifies  the  solar  radiation  reaching  the  earth,  due  to  both 
scattering  and  cibsorption.  Characterization  of  optical 
properties  of  atmospheric  aerosols  with  the  proposed 
nephelcmeter  will  provide  valuable  information  needed  to  develop 
models  describing  the  radiative  energy  exchange  in  the  earth- 
atmosphere  systan,  and  the  effects  which  the  addition  of 
aerosols  has  on  this  exchange.  Quantitative  analysis  of 
atmospheric  scattering  phenomena  are  often  limited  by  the  fact 
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that  little  information  is  available  about  refactive  index 
values.  Differential  scattering  data  can  be  used  to  determine 
the  conplex  index  of  refraction  of  airborne  particles. 

(5)  Optical  Particle  Sizing; 

Optical  particle  sizing  instruments  developed  to  date  give 
erroneous  results  for  nonspherical  particles  application.  Based 
on  the  knowledge  obtained  from  measuranents  with  the  proposed 
nephelometer ,  development  of  a  new  instrument,  which  can  provide 
an  in-situ  measurement  of  particle  size  of  nonspherical 
particles,  can  be  directed. 

5.2  Potential  Use  by  the  Federal  Government 

Potential  uses  of  the  multi-wavelength,  multichannel  nephelometer 
system  developed  through  the  proposed  project  are  expected  in  the  field 
of; 

(1)  Performance  Analysis  of  Optical  and  Infrared  Systems; 

Performance  of  optical/infrared  systems  such  as;  thermal 

imaging  system,  night  vision  devices,  FLIR  system,  guidance/seeker, 
laser  rangefinder,  laser  designator,  optical  communication,  and 
directed  energy  weapon  system  is  greatly  affected  by  atmospheric 
conditions  under  which  those  systems  have  to  be  operated. 
Characterization  of  optical  properties  of  atmospheric  aerosols  with 
the  proposed  nephelometer  will  enhance  one's  ability  to  predict  EO/IR 
systems  performance  under  different  atmospheric  conditions. 

(2)  Optical  Countermeasures ; 


Characterization  of  optical  properties  of  aerosol  particles  will 
provide  an  understanding  of  principles  governing  the  optical  merit 
of  obscurants.  Effectiveness  of  optical  countermeasures  such  as 


obscuTcuit  smokes  can  be  analyzed  based  on  the  phase  function  data 
measured  with  the  proposed  nephelometer. 


(3)  Chenical/Bioloq ical  Agent  Detection; 

Correlaticai  of  particle  optical  properties  measured  with  the 
proposed  nephelometer  with  known  chemical  and/or  biological 
constituents  will  greatly  advance  the  development  of  CTJ  agent 
detection  system. 

(4)  Atmospheric  Modeling  of  Optical/Infrared  Radiation  Transmissions 
Characterization  of  optical  properties  of  atmospheric  aerosols 

can  provide  a  data  base  for  an  attempt  to  model  and  predict  their 
effects  on  the  propagation  of  optical  and  infrared  radiation.  These 
effects  include  transmission  loss,  contrast  degradation,  and  path 
radiance.  Optical  paths  through  the  atmosphere  include  horizontal 
and/or  slcint  paths. 

(5)  Remote  Sensing 

Knowledge  on  the  optical  properties  of  aerosol  particles  will 
greatly  enhance  ones  ability  to  interpret  the  experimental  data 
measured  with  remote-sensing  instruments  such  as  multi-wavelength 
solar  radiometer,  monostatic  and  bistastic  lidars,  and  imaging 
systems. 
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6,0  Ccxiclusions 


Design  parameters  for  a  multi-wavelength,  multichannel  nephelaneter 
system,  vAiich  is  capable  of  measuring  differential  scattering  patterns  of 
individual  aerosol  paticles,  have  been  optimized.  System  design  includes  the 
particle  saitpling  inlet,  the  optical  system,  and  the  data  acquisition  system. 
The  new  system  provides  means  for  making  in-si tu,  real-time,  multi-wavelength 
measurements  of  the  optical  properties  of  individual  aerosol  particles  of 
different  shape  and  composition.  ^stem  uses  0.6328,  3.39,  and  10.6  pm 
laser  wavelengths,  v*iich  are  in  the  visible  and  infrared  atmospheric  window 
regions.  The  optical  system  measures  multi-wavelength  differential  scattering 
patterns  of  a  single  particle  in  three  different  scattering  planes  to  account 
for  effects  of  particle  orientation  on  scattering  by  nonspherical  particles. 
The  particle  sampling  inlet/systan  sanples  particles  with  minimal  dependence  on 
the  micrcmeteorological  conditions  of  the  ambient  atmosphere  and  directs  them 
in  a  train  of  single  particles  into  the  scattering  volume.  The  ccmputer-based 
data  acquisition  system  can  acquire  the  measured  differential  scattering  data 
and  store  than  simultaneously  into  a  mass  storage  device  without  any  data  loss. 

The  Phase  II  research  and  development  effort  involves  fabrication, 
calibration  and  performance  testing  of  the  multi-wavelength,  multichannel 
nephelcroeter  based  on  the  specific  design  of  the  system  resulting  from  the 
Phase  I  research  work.  This  new  nephelometer  system,  if  it  is  developed 
successfully  during  the  Phase  II  effort,  will  be  able  to  provide  information 
which  can  enhance  one's  knowledge  concerning  the  optical  properties  of  airborne 
particles.  If  this  knowledge  is  obtained,  then  the  ability  to  predict 
performance  of  optical  and  infrared  systans  through  the  atmosphere  is  enhanced. 
The  uniqueness  of  this  approach  lies  in  the  inherent  ability  of  the  system  to 
produce  data  heretofore  unattainable. 


The  benefits  resulting  from  the  successful  development  of  such  an 
instrument  are  seen  in  numerous  areas  of  military  and  civilian  applications. 
In  general,  researchers  and  system  designers  will  have  an  instrument  which  will 
provide  more  information  'than  presently  available  with  respect  to 
optical/physical  properties  of  individual  atmospheric  and/or  artificial 
particles.  Upon  completion  of  development  and  performance  testing  of  a 
prototype  of  the  proposed  nephelometer  systau  during  Phase  II  effort, 
comnercial  ajplicaton  of  such  a  systan  will  be  initiated  during  Phase  III 
effort. 


7.0  Participating  Professional  Personnel 


FAME  Associates,  Inc.  members  who  participated  in  Fhase  I  research  effort 
are  listed  below: 

James  B.  Wedding,  Ph.D.,  Principal  Investigator 
Young  J.  Kim,  Ph.D.,  Senior  Research  Engineer 
Michael  A.  Weigand,  Director,  Aerosol  Science  Laboratory 
None  of  the  project  staff  earned  advanced  degrees  during  the  course  of  the 
project. 
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